NASA Contractor Report 194965 
ICASE Report No. 94-71 



ICASE 


MODELING JETS IN CROSS FLOW 

(NASA-CR-19A965) MODELING JETS IN 
CROSS FLOW Final Report (ICASE) 

37 p 



G3/34 


A.O. Demuren 


Contract NAS 1 - 1 9480 
August 1994 


Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, VA 23681-0001 


N95-1 1930 
Unc las 
0022778 


SRA Operated by Universities Space Research Association 



MODELING JETS IN CROSS FLOW 


A.O. Demuren 1 

Department of Mechanical Engineering 
Old Dominion University 
Norfolk, VA 23529 


ASTRACT 

Various approaches to the modeling of jets in cross flow are reviewed. These are grouped 
into four broad classes, namely: empirical models, integral models, perturbation models, and 
numerical models. Empirical models depend largely on the correlation of experimental data and 
are mostly useful for first-order estimates of global properties such as jet trajectory and velocity 
and temperature decay rates. Integral models are based on some ordinary-differential form of the 
conservation laws, but require substantial empirical calibration. They allow more details of the 
flow field to be obtained; simpler versions have to assume similarity of velocity and temperature 
profiles, but more sophisticated ones can actually calculate these profiles. Perturbation models 
require little empirical input, but the need for small parameters to ensure convergent expansions 
limits their application to either the near-field or the far-field. Therefore, they are mostly useful 
for the study of flow physics. Numerical models are based on conservation laws in partial- 
differential form. They require little empirical input and have the widest range of applicability. 
They also require the most computational resources. Although many qualitative and quantitative 
features of jets in cross flow have been predicted with numerical models, many issues affecting 
accuracy such as grid resolution and turbulence model are not completely resolved. 
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Introduction 


In the earliest studies of jets in cross flow, empirical models were developed to correlate 
experimental data obtained under various idealized conditions. Such models are reviewed in 
detail by Abramovich (1963), Rajaratnam (1976) and Schetz (1980). They mostly give the 
jet trajectory and center-line decay rates. The earliest approach based on the actual solution 
of conservation equations belongs to the class of integral models. These are derived either 
by the application of conservation principles to a finite control volume or by the use of profile 
assumptions to simplify the partial differential equations which describe conservation laws. These 
models offer more flexibility than empirical models. Jet trajectories, decay rates, growth rates, 
and even cross-sectional shape have been predicted. However, empirical input is usually required 
in the form of entrainment rates or drag coefficient. Further, it may be difficult to prescribe cross- 
sectional profiles in complex situations. Numerical models attempt to solve some form of the 
full partial differential equations, which represent conservation laws, by using a finite-difference, 
finite-volume or finite-element method. Little empirical input is required, hence they have the 
potential for the widest range of applicability. However, there may be problems with inadequate 
grid resolution, imprecise boundary conditions and deficiencies in the turbulence model used for 
closure of the mean flow equations. Recently, several models based on perturbation methods 
have been proposed. These are mostly of scientific interest, since drastic assumptions such as 
inviscid flow, negligible jet distortion, small deflection, etc., may be required for the perturbation 
analysis. Thus, they are used mainly to study the flow physics in limited regimes, either in the 
near-field or in the far-field. 

In practical engineering applications, jets in cross flow are found in both confined and 
unconfined environments. Examples of confined jets in cross flow include: 1) Vertical and Short 
Take-Off and Landing (V/STOL) aircraft in transition from hover to forward flight, in which case, 
the jets from its engines impinge on the ground surface; 2) Internal cooling of turbine blades 
by air jets impinging on the leading edge, and; 3) Dilution air jets in combustion chambers of 
gas-turbine engines, where the jets are injected radially into the chamber, through discrete holes 
along its circumference, in order to stabilize the combustion process near the head, and to dilute 
the hot combustion products near the end. 



Practical examples of jets in unconfined or semi-infinite cross flow are more numerous. 
These include: 1) Flow situations resulting from the action of cross winds on effluents from 
cooling towers, chimney stacks, or flames from petrochemical plants; 2) Discharge of sewage 
or waste heat into rivers or oceans; 3) Film cooling of turbine blades; 4) The use of air curtains 
to prevent cold air from entering open spaces in industrial buildings, and; 5) Thermal plumes 
rising into cross winds in the atmosphere. 

The configuration of a jet in cross flow is illustrated in Fig. 1. The axis of the jet is usually 
defined as the locus of the maximum velocity or total pressure. The jet trajectory is referred 
to this line, as opposed to the center-line of the jet, which is mid-way between the inner and 
outer boundaries of the jet, usually determined from flow visualization. The main parameter 
which characterizes a jet in cross flow is the jet-to-cross-flow velocity ratio, R (= Uj/U 0 ), or the 
momentum flux ratio J (= ^-R 2 ). In confined jets, the normalized wall distance H/D may also 
be important if it is not very large. In multiple jets, the normalized spacing S/D will be a factor. 

As shown in Fig. 1, the jet in a cross flow has three main regions: the potential core zone 
(I); the zone (II) of maximum deflection; and the far-field zone (HI). The potential core, in the 
central part of zone I remains relatively unaffected by the cross flow though its length is reduced 
in comparison to that of a jet in stagnant surroundings. Thus, for a turbulent jet it reduces from 
(~ 6D) to 6.2D e~ 3 3/R [Fan (1967), Pratte and Baines (1967)] or 6.4/(l+4.6/R) [Kamotani and 
Greber (1972)]. The two relations deviate at low R, where the potential core length is strongly 
influenced by actual exit flow conditions. In zone II, the jet experiences the most deflection. The 
pressure gradient across the jet is maximum as well as the entrainment rate. This is the most 
difficult zone to analyze accurately. In the far-field zone ID, the jet axis approaches the crossflow 
direction asymptotically, and the flow field is nearly self-similar. All models can predict this 
region fairly well, given the correct boundary conditions at the end of zone II. 

Four broad classes of models, namely, empirical, integral, perturbation and numerical, are 
now described. Emphasis will be given to the last two, where most of the recent advances have 
been. The first two methods were reviewed extensively in monographs by Abramovich (1963), 
Rajaratnam (1976) and Schetz (1980), and in a review article by Demuren (1985a). 
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Empirical Models 

Empirical models present the simplest means of predicting global properties of jets in cross 
flow. They depend largely on the correlation of experimental data, and the accuracy of the 
predictions may depend on the closeness of the conditions of the particular problem of interest 
to those in the data base used for the correlation. Due to their low cost and ease of use, empirical 
models are most useful for first-order estimates and as qualitative checks for results produced 
by other methods. 

The most common parameter given by empirical models is the jet trajectory. For a single 
circular turbulent jet injected normally into a cross flow, the trajectory has the form: 



where, in the range of J between 2 and 2,000, a has a value between 0.7 and 1.3, b has a value 
between 0.36 and 0.52 and c takes a value between 0.28 and 0.40, depending on experimental 
conditions. The values a=0.85, b=0.47 and c=0.36 appear to be a good compromise for the 
intermediate range of J. This equation should also be valid for confined jets, up to the point 
of contact, and for multiple jets with medium to large spacing ratios. Equation (1) with b = 
0.36, and c = 0.28 also gives the physical boundaries of jets in cross flow [Pratte and Baines 
(1967)]; with a = 1.35 and 2.63 for the inner and outer boundaries, respectively, and a = 2.05 
for the center line. For plane jets in confined cross flow Kamotani and Greber (1974) found 
that Eq. (1) can be used with a = 2.0 (l-e _H/D ), b = 0.28 and c = 0.50. Equations for other 
parameters such as entrainment rates, velocity profiles, temperature trajectories, etc., can be 
found in Demuren (1985a). 

Integral Models 

Integral models are the first elaborate calculation procedures applied to predict the behavior 
of jets in cross flow. In these models, integral equations are derived either by considering 
the balance of forces and momentum changes over an elementary control volume of the jet, 
or by integrating in two spatial directions, the three-dimensional, partial differential equations 
governing the jet flow. In either case, a set of ordinary differential equations is obtained which 



can be solved analytically or numerically. Empirical input is required to prescribe pressure drag, 
entrainment rates and spread rates. The former approach is easier to understand and to implement 
and is therefore more popular. On the other hand, the latter approach involves more extensive 
mathematical manipulation, but it is more transparent in the assumptions made and affords more 
flexibility in dealing with complex boundary conditions and trajectories. 

Integral models flourished between the late 1960’s, when more flexibility was required 
than could be obtained with empirical models and the early 1980’s, when the rapid growth in 
computer hardware and software made elaborate numerical computations of three-dimensional 
flows feasible. Many of the earlier model are reviewed by Rajaratnam (1976). In these models, 
there was an assumption of the constancy of the momentum in either the initial jet direction, the 
cross flow direction or the axial direction, and the jet was bent over by a prescribed pressure 
drag force, or entrainment of ambient fluid. None of these models could predict correctly the 
jet trajectory over a range of R [Demuren (1985a)]. Thus, they offer no advantage over much 
simpler empirical models. 

More refined integral models consider effects of both the pressure drag and the entrainment 
of cross flow ambient fluid on the jet. A typical model is that proposed by Fan (1967) for buoyant 
jets in cross flow. The cross section of the jet was assumed circular with radius y/2 b, and the ex- 
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cess velocity profile was assumed to be Gaussian, i.e., V - U D cos 6 = (V max - U 0 cos 8)e~ T> ' . 
The resulting set of ordinary differential equations can be written as: 

Continuity 

^(pAV)=Cp 0 U, 


x-momentum 


y-momentum 


— (pAA v V 2 COS0) = C p 0 U e U 0 + 0.5 CdAz p 0 U 2 sin 3 6 


— (/5AA v V 2 sin0) = -A (p - p 0 )g + 0.5 CqAz p 0 U 2 sin 2 8 cos 0 


Scalar (Temperature/Concentration) <f> 


dC 


(p A A <j>) = 0 
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where A is the cross-sectional area of the jet, C its circumference, p is the density of the jet 
fluid, g the acceleration due to gravity, and A v and A^ are respectively, momentum or scalar flux 
coefficients which depend on the assumed velocity and scalar profiles. U e is the entrainment 
velocity. Fan proposed that it should be proportional to the velocity vector difference between 
the jet and the cross flow, but this was found to be unreliable. Abraham (1971) proposed an 
entrainment model with two parts as: 

Ue = E mom (V max U 0 cos 0) -f- Ef^Uo sin 6 cos 0 

where the coefficients E m0 n and E t h have the values 0.057 and 0.50, respectively. The first part of 
Eq. (6) represented the entrainment of a momentum jet in a nearly stagnant ambient fluid and the 
second part the entrainment into (momentum-free) thermals under similar conditions. Cos# was 
introduced artificially into the second part to prevent it from contributing to entrainment when the 
jet was nearly perpendicular to the cross flow. With this entrainment model, the drag coefficient 
Cd was given a value of 0.3. Equations (2) through (6) were then integrated numerically. With 
this model, Abraham (1971) was able to obtain quite good agreement with experimental data 
of jet trajectory and axial concentration decay. Typical results are shown in Fig. 2. For a 
non-buoyant jet, E t h = 0 and p = p 0 so that the buoyancy term in Eq. (4) is also zero. 

Equations (2) through (6) may also be applied to predict plane jets in cross flow by 
substituting the appropriate expressions for the area A and the circumference C. The entrainment 
coefficients E m0 m and E t h and the drag coefficient Cd must then be calibrated with plane jet data. 

In order to extend the range of applicability of integral models to include more flow physics 
and be able to deal with more practical situations, such as multiple jets in varied arrangements, 
highly non-uniform cross flow, etc., more elaborate models have been proposed by Campbell 
and Schetz (1973), Isaac and Schetz (1982), Makihata and Miyai (1983), amongst others. All 
these models were derived based on the control volume approach and are applicable mainly to 
jets with plane trajectories. Hirst (1972) and Schatzmann (1979) developed models based on the 
integration in the cross plane of the jet of the three-dimensional partial differential equations by 
making the assumption of axi-symmetry and that profiles of the excess velocity are Gaussian. 
These latter models could be applied to situations with three-dimensional jet trajectories. 
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Although integral models allow economical prediction of several flow properties in com- 
parison to full-blown numerical models, they have been criticized for the need to assume the 
shape of the jet cross-section and profile functions, some of which may not be realistic for the 
whole evolution of the jet in cross flow, especially in the zone of maximum deflection. How- 
ever, it appears that in spite of the apparent oversimplification, integral models can be made to 
perform well in some cases with proper calibration. Adler and Baron (1979) proposed a quasi- 
three-dimensional model which did not assume the cross-sectional shape or similarity profiles 
for the velocity, but these were computed along with other flow variables. The characteristic 
kidney-shaped cross-section of the jet was computed successfully by considering the evolution of 
vortices distributed along the boundaries of the jet in a Lagrangian manner, and the cross-section 
was allowed to grow at a rate which was an average between the growth rates of free jets and 
vortex pairs. Similarly, velocity profiles were allowed to change in the zone of maximum deflec- 
tion, culminating in self-similar profiles only in the far-field zone (HI). The model gives quite 
good prediction of the three-dimensional flow fields of jets in cross flow studied experimentally 
by Kamotani and Greber (1972). 

Perturbation Models 

If in the jet in cross flow problem a small parameter is defined, perturbation methods can 
be used to solve the governing equations. Most applications of perturbation models have been 
to study the flow of strong jets in a weak cross flow. In the initial stage, the flow can be 
considered to be a small perturbation from that of a free jet in stagnant surroundings, and the 
jet stiffness A (= 1/R) can be used as the small parameter. This places a severe restriction on 
the range of applicability of such models. However, they have the advantage of not being too 
dependent on empirical calibration as integral methods are, and they are computationally much 
cheaper than numerical methods. The goal is to predict the main features of jets at high R ( 
>10), including trajectory, cross-sectional shape, velocity field, vorticity field, mixing, etc., with 
minimal empirical input. It was believed [Needham et al . (1988), (1990)] that the jet distortion 
and deflection could be obtained by inviscid analyses based on the evolution of vortex filaments 
around the jet as it exits from a circular pipe or orifice. This approach was based on the 
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original work of Chen (1942) which was also the basis of fairly successful computations with 
the integral model of Adler and Baron (1979). 

Chen’s model approximates the near-field as two regions of irrotational flow, the jet flow 
and the external cross flow, separated by a vortex sheet. The three-dimensional vortex sheet 
can then be approximated by a two-dimensional vortex sheet, which originates from the pipe 
or orifice exit and evolves in time in the axial direction. Needham et al . (1988), (1990) 
applied a three-dimensional model, with perturbation expansions for the potential flow within 
and outside of the jet. The distortion of the jet could be predicted reasonably well, but contrary 
to earlier studies, no jet deflection was obtained, if the jet issued normally into the cross flow. 
Surprisingly, with a component of the cross flow in the direction of the jet, some deflection 
was obtained. This discrepancy was explained by Coelho and Hunt (1989) who showed that 
the two-dimensional time-evolving vortex sheet model was a poor approximation for the fully 
three-dimensional vortex sheet model. 

The two-dimensional vortex sheet equation can be written as 


W + 5J ((/ ' 7) "° 


where 7 is the vortex strength and U s is the average speed of the flow across the layer, using 
the nomenclature of Fig. 3(a). If 7 (9, t), U s ( 9 , t) and R(0, t) are approximated as Taylor series 
expansions with respect to t, derivatives of 7, U s and R of any order with respect to t, at t=0, 
can be evaluated. This gives the shape of the vortex sheet or the jet boundary as 


R(9,t) = R 0 - 


n u 2 i 

--f cos 20 

t 2 + ’ 

2 Ro 



1 U 3 

--^-{3 cos 3# — cost?} 
6 Rl 


t 3 + 0(t 4 ) 


2y 


If it is assumed that elements of the vortex-sheet travel at half-speed, such that t = — , and Uj 

U i 


and Ro are used for normalization, then Eq. (8) becomes: 


R(9,y) = 1 - [2*/ 2 cos 29] A 2 - 


— j/ 3 { 3 cos 3 9 — cos 9 } 

O 


a 3 + o(a 4 ) 


However, in the presence of the cross flow, the flow in the jet pipe is distorted [Andreopoulos 
(1983)], so that Uj is not uniform, and Eq. (9) may not be a good approximation for Eq. (8). 
The actual non-uniformity of the jet exit flow can be calculated using a fully three-dimensional 
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vortex-sheet model. In this case, the longitudinal and transverse components of vorticity must 
be considered. As shown in Fig. 3(b), these can be approximated by the vertical and azimuthal 
components 7 y and 7 S , respectively. The general expression for the strength of the vortex sheet 
can be written in vector form as 


^7+ (UvV) 7 = (7V)U V --?(V-U r ) 


where 7 = { 7 3 , 7 ;,},U V = {U s ,U y } and V 


Ids’dt/J- 


The vertical component of Eq. (10) is 




dUy 

ds 


'Is 


which contains a source term, in contrast to Eq. (7). This source term expresses the rate at which 
fluid elements rotate as they travel up the vortex sheet. Thus, the vertical vortex strength may 
be strongly influenced by the azimuthal vortex strength and variations in the azimuthal velocity. 
The solution of the three-dimensional vortex-sheet problem, in terms of the potential flow inside 
the jet, and the external potential flow gives the shape of the vortex-sheet to third order as 


R{0,y) = l-A 2 


y 2 - 2 C 2 y - AnH^n)e- 9 ^ - 1 

+0(A 4 ) 


71=1 


cos 29 


where C 2 is a constant and A n are coefficients given by the boundary conditions. J 3 ( ) are 
Bessel functions of third order, and cr n are zeros of J 2 ( ). Solutions for the velocity field (in 
terms of tlje velocity potential) and the pressure field are also given in terms of Bessel functions. 
Comparison of Eqs. (9) and (12) shows that the O (A 3 ) term in the former, which produces the 
deviation from symmetry and thus the jet deflection is absent in the latter. Therefore, Coelho and 
Hunt (1989) concluded that a three-dimensional inviscid vortex-sheet model could not produce 
jet deflection. Although, the two-dimensional model of Chen (1942) could produce a deflection, 
this was an incorrect approximation of the three-dimensional flow. However, by introducing a 
vertical component to the cross flow, the symmetry in Eq. (12) may be broken and a second 
small parameter is introduced into the perturbation expansion, as in the works of Needham et 
al. (1988), (1990). Then, jet deflection would occur. 
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Coelho and Hunt (1989) postulated that viscous or turbulent entrainment was necessary for 
jet deflection. They proposed an entraining vortex-sheet model [see Fig. 3 (c) for nomenclature]. 
The entrainment velocity from the external cross flow is given by 

K, = e( 7j 2 +7 , 2 )' /2 + 0(e 2 ) 


where e is the entrainment coefficient, which must be prescribed empirically, and it now becomes 
a second small parameter for the perturbation expansion. The mixing layer within the vortex- 
sheet entrains fluid from both the jet and the external cross flow. To the leading order, the 
entrainment rates are assumed proportional, so tjiat 

Vje = e c{j 2 s + 7 ^) 1/2 + 0(e 2 ) 


where c is a constant of O (1). The equation for the conservation of mass within the mixing 
layer is 


|;[<W 


R ’ )] + ^ 


Uy 



e (ls + ly) 


1/2 




+ C 


Rj + 



+ 0(e 2 ) 


Equation (15) must be solved along with the potential flow equations for the jet flow and the 
external flow. The solution yields for the mean radius 

R(0, y) = 1 + ye + [(1 + c)y 2 cosO ] eA — [Z(y)cos26}\ 2 

+0(e 2 , A 3 , eA 2 ) 

OC 

where Z(y) = y 2 — 2 C 22 / — A n J$(<r n ) [e~ anV — l]. Comparison of eqs. (12) and (16) shows 

n= 1 

that the deviation from symmetry is now of O(eA), so that jet deflection would occur, as one 
would normally expect. 

Higuera and Martinez (1993) have proposed a mixed perturbation/numerical model which 
does not use the vortex-sheet concept but solves the parabolized Navier-Stokes (PNS) equations 
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in the distortion region of jets in weak cross flow (R>15). The model is applicable to laminar flow 
or a turbulent flow in which the assumption of a constant eddy viscosity would be appropriate. 
The weak cross flow is necessary so that there is only mild curvature in the distortion zone II, 
enabling the governing equations to be parabolized. Furthermore, there should be little deviation 
of the jet flow in the development zone I from that of a free jet, so that the flow field at the end 
of zone I and the beginning of zone II can be prescribed from Landau’s self-similar solution for a 
point source of momentum [Batchelor (1967)]. A perturbation method is used to solve the PNS 
equations for small y, with y as the small parameter. For intermediate values of y, a parabolic 
numerical method is used. However, computations must be stopped once the distortion becomes 
too large for the assumptions of negligible axial diffusion and pressure gradient to be valid. 

For small y = [0(4y w )], where y w = [(\/3/16)RD], the deflection of the jet, ft, will be small. 
The continuity and momentum equations, in dimensionless variables, can then be written for the 
axial velocity component, v and cross stream velocity vector V = [«,w] as 


-» -> dv 

V ■ V + — = 0 

dy 


dv (- al( 3P 0 2 v\ 

v^+{V-V)v = V v + fl + 


dv , _ ^ 

V— + (V • VJ V = -VP + V 2 V + ft' 


, d 2 V 

dy 2 


where V = { 57 } and V 2 = The requirement that y w should be beyond the 

development zone I indicates how large R must be for the analysis to be valid. For example for 
turbulent flow, with a development length «6D, R > 15. Therefore ft — > 0, and terms in ft 2 can 
be neglected. If the pressure gradient term is also eliminated by combining the divergence of 
Eq. (19), with the continuity equation, the governing equations become 


£ 2 + (*.*)—*», 


dy 


dfl / - - \ dv du dv dw dv 2 

v-v )n = n— + + v 2 o 

ay dy oz dy dx 


— + 6 
dy + l 


dw_du 
dx dz 
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where SI is the vertical component of vorticity. Equations (20) - ( 22 ) are parabolic in y, so they 
can be solved by marching in zone II. 


The initial conditions are derived from Landau’s self-similar profiles as 

( yu , yw, yv, y 2 fi) — »• ( V rs cos 6, V r3 sin0, v 3 , 0) as y — ► 0 

4t/(1 — y 2 ) 8 r 

9“ , v s = 9 , 77 = - 

(1 + 7 ] 2 ) (1 + T} 2 ) V 

-> 00 are 

v = w = S} = 0 ; u — ► 1 


where. 


The boundary conditions as x 


the near-field solution has the form 


(v,v r ,v e ,ySl) = -(v ^,0,0) -I- (V 0 ) ,t4 0 ) ,u ( 


( 0 ) 

6 



K (1) , -l- 

The functions on the right hand side of Eq. (25) depend only on 77 and 6, whereas Eqs. (20) - (22) 
are functions of 77 , 6 and y. Hence, by substituting Eq. (25) into Eqs. (20) - (22) and collecting 
terms of like order in y, solutions of different order can be found. Of course, terms of order 
(-1) will reproduce the initial conditions. This approach is really quite restrictive. The several 
requirements of large R, small y, constant turbulent eddy viscosity and low Reynolds number 
exclude it from consideration as a realistic tool for practical computations of jets in cross flow. 

In general, perturbation models are not yet sufficiently mature to become more than curious 
tools of analysis. A redeeming factor is that it is especially in those high R flows, for which they 
are valid, that most numerical models are least accurate. In these flows, there are substantial 
regions with high shear and rates of strain, in which standard discretization schemes and 
turbulence models may become inadequate. 
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Numerical Models 


Numerical models have the most potential for wide generality and can, in principle, be 
applied to the whole range of jet in cross flow situations, confined or unconfined, low medium 
or high R, single or multiple jets, impinging on a wall or on other jets, swirling, homogeneous 
or heterogeneous cross flow, compressible or incompressible, etc. The analysis starts from the 
general conservation laws stated in partial differential equation form, which are the Navier- 
Stokes equations for the velocity field, and corresponding energy or species equations for the 
temperature or concentration fields, respectively. These equations, which describe unsteady, 
three-dimensional flow cannot be solved directly in practical applications for turbulent flows. 
In incompressible fluid flow, time-averaged forms, and in compressible fluid flow, density- 
weighted, time-averaged (or Favre-averaged) forms of the equations are solved. The process of 
time-averaging introduces a closure problem due to non-linear correlation between fluctuating 
velocity and/or temperature/concentration fields. Turbulence models are required to determine 
these correlations, thereby affecting closure of the system of equations. Most numerical models 
applied to the jet ia cross flow problem use the eddy viscosity concept. In its simplest form, the 
turbulent eddy viscosity is prescribed as a constant, whereas more sophisticated models solve 
partial differential equations for turbulent quantities, from which the eddy viscosity distribution 
can then be obtained. Experimental studies by Andreopoulos and Rodi (1984) show that there 
are significant regions of the jet cross flow interactions in which the eddy viscosity concept is 
invalid. Demuren (1992) proposed a numerical model in which the eddy viscosity concept is not 
invoked but partial differential equations are solved to determine distributions of the turbulent 
correlations directly. In most numerical models, the computational domain encompasses the 
whole region in which the influence of the jet is felt, or if necessary the whole field of the jet 
and cross flow. No assumptions are required as to the evolution of the jet within the flow domain, 
but this is obtained as a result of the computations. It is only necessary to prescribe boundary 
conditions at the chosen computational boundaries. The two major issues in the application of 
numerical models to jets in cross flow are the accuracy of the basic numerical method and the 
accuracy of the turbulence model. 
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The time-averaged, three-dimensional, steady-state mean flow equations can be written in 
Cartesian tensor notation as 

continuity 


d_ 

dx\ 


(pUt) = 0 


momentum 


h (pUiUl) = -t, p+ k 


■ (<>U, (III, 


d : 


scalar 




—p u\d + 


p d<f> 

a dx[ 


with i = 1, 2, 3 and / = 1, 2, 3 representing properties in the lateral, vertical and longitudinal 
directions, respectively. The equations are expanded with Einstein’s summation rule for repeated 
indices, x, are the Cartesian coordinates and C/, the Cartesian velocity components. 4> may 
represent any scalar such as the temperature or species concentration, —pupal and — pu\d 
represent the Reynolds stresses and the turbulent scalar fluxes, respectively. Distributions of 
these quantities are obtained from the turbulence model. Also, p is the molecular viscosity and 
a the corresponding Prandtl of Schmidt number. is the source term for the temperature or 
concentration equation. 

The task of the turbulence model is to provide distributions for the Reynolds stresses and 
the scalar fluxes so that the mean flow Eqs. (26) to (28) can be closed. In the Boussinesq eddy 
viscosity concept, the Reynolds stresses are calculated from 


-pu^i = p t 



pkS tt 


The corresponding eddy diffusivity concept gives 


— p U/t) 


Pt 

crj, dxi 


where p t is the turbulent eddy viscosity, a ^ is the turbulent Prandtl or Schmidt number, k is the 
turbulent kinetic energy (per unit mass) and 8u is the Kronecker delta which is equal to unity 
when i = l, and zero otherwise. 
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The most common method for calculating the distribution of pt is through the k-e turbulence 
model [Launder and Spalding (1974)]. This gives 

k 2 

(H = p — 

e 


The distributions of k and e are then obtained from solution of transport equations which can 
be written in Cartesian tensor form as 


m + ? Pt - 1,1 




dxi 


d ( pt de \ 
dxi \<T t dx\ ) 


n e e 

” 1 " C €lP _ 


where e is the rate of dissipation of k, and P* is the rate of production of k through the interaction 
of the Reynolds stresses with the mean flow. It is given by 

. dUi 


Pk = -%«i 


dx. 


The empirical coefficients which appear in Eqs. (29) — (33) are given the standard values Cp 
= 0.09, c £l = 1.44, c t2 = 1.92, <r ^ = 0.9, crjj = 1.0 and <r t = 1.3. 

Simpler eddy viscosity relations have been utilized with reasonable success in some studies. 
Chien and Schetz (1975) prescribed a constant value for p t , proportional to the jet velocity 
excess and the jet diameter. However, in a subsequent study Oh and Schetz (1990), calculated 
p t from a relation which takes into consideration the complex shape of the jet and the relative 
magnitude of the axial turbulence intensity to the velocity excess. Thus 


pt — 0.037 / p fej/2 Af/ C 


where bi# is the characteristic half width of the jet, A U c the centerline velocity excess, and / 
( = u' 2 /AUc) takes a value of 2 in the potential core and {1 + exp[— 1.134(C - ( 0 )]} in the 
main jet region. For computations at high R, Sykes etal. (1986) proposed a one-equation model 
which solves the k equation such as Eq. (32), but calculates the Reynolds stresses from 
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where A is a length scale given by 


A = 0.088D + 0.0088r 


D is the jet diameter and r is the distance from the center. In spite of the rather crude 
length scale assumption, computed results, of the mean flow (for R=2) agreed reasonably 
with experimental data. However, turbulent kinetic energy levels were grossly overpredicted, 
especially in the wake region. 

Demuren (1992), (1994) and Alvarez and Jones (1993) have used various Reynolds stress 
models (RSM) to investigate the effect of the turbulence model on computations of jets in 
cross flow. Model computations in which the eddy viscosity concept is not invoked but partial 
differential equations are solved for the Reynolds stresses are compared to those using the k-e 
model and to experimental data. The Reynolds stress equations can be written in Cartesian 
tensor notation as 


d_ 

dx} 


(U\ u t u } ) — D(j -(- P{j -T Ttij f ij 


where D,y is the turbulent diffusion, P tJ is the production, 7r,y is the pressure-strain correlation and 
e,y the dissipation rate. The production term P,j = —uiui ^ — ujuj and the dissipation 
rate is assumed to be locally isotropic so that =2/3 S tJ e. D,y and 7r,y contain higher-order 
correlations, and so must be approximated for closure at this level. In Demuren (1992), these 
terms are modeled after proposals of Daly and Harlow (1970) (denoted by DH) and Launder, 
Reece and Rodi (1975) (denoted by LRR), respectively. In Demuren (1994) and Alvarez and 
Jones (1993), additional models for D,y and 7r,y are considered, including those proposed by 
Mellor and Herring (1973) (denoted by MH) and Speziale, Sarkar and Gatski (1991) (denoted 
by SSG), respectively. The latter combination of models was found to give the best overall 
predictions of developed turbulent plane channel flow [Demuren and Sarkar (1993)]. The DH 
and MH diffusion models can be written, respectively, as 


Dij = 


C$2 J~ 

dx k 


k 2 


Oil ill j 


+ 


e \ dxk dxj dxi J\ 
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with C s i = 0.22 and C s 2 = 0.072. The pressure-strain models can be written in the general form 


ij — a o e by -j- cx\ c^bj k b k j ^ll^y^ + 02 k Sy 

“b ^3 F*k by "T k{a 4 ^bj k Sj k -F bj k Sj k ^^ijbkl 
+05 (b; k Wj k -|- bjkWjk)} 


where by = | 5y j is the Reynolds stress anisotropy tensor, Sy = \ j is 

the rate of strain tensor, W;j = is the rotational tensor, and II = bj k b k / is the 

second invariant of by. For the LRR model, a G = — 3 + f w , ari = 03 = 0, «2 = 0.8, a 4 = 1.745 
and 05 = 1.309 — 0.24 f w . For the SSG model, o Q = —3.4, oj = 4.2, 02 = 0.8 — 1.3 II 1 / 2 , 
03 = — 1.8, o 4 = 1.25 and 05 = 0.40. f w is a wall proximity function which has a value of 
unity near a wall and zero in a turbulent flow free from walls. The correct rate of decay away 
from walls is a subject of controversy [Demuren and Rodi (1984)]. It is also difficult to specify 
in complex flows with curved walls or multiple walls. The absence of such a term makes the 
SSG pressure-strain model rather attractive for application to complex flows. 

A turbulence modeling approach which is intermediate between the k-e model and the full 
Reynolds stress model was utilized by Baker et al . (1987) to calculate the near field of jets in 
cross flow at high R. This is the so-called algebraic Reynolds stress model, which is derived 
by dropping the convection and diffusion terms in Eq. (38). Thus, implicit algebraic equations 
are obtained which can be solved simultaneously for the Reynolds stresses. Baker et al . (1987) 
used further simplifications of these equations to obtain explicit expressions for the Reynolds 
stresses. These expressions contain k and e as unknowns so that Eqs. (32) and (33) must still 
be solved before the Reynolds stresses can be calculated. This approach falls within the general 
class of non-linear k-e models reviewed by Speziale (1991). 

The boundary conditions applied to the equations depend on the particular problem. The 
various types of boundaries which may exist in these flow situations are inflow, outflow, wall, 
symmetry planes and the free stream. At inflow boundaries, the values of the dependent variables 
are prescribed or deduced from experimental data. At outflow planes, boundary conditions such 
as zero traction force or zero normal gradient are usual prescribed. It is normal to prescribe 
no-slip conditions along walls, but in order to bridge the flow between the fully turbulent region 
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and the viscous sublayer near the wall, the wall-function method [Launder and Spalding (1974)], 
which is based on the assumption of local equilibrium, is usually employed to prescribe variable 
values along the first set of grid nodes nearest to the wall. Kim and Benson (1992) did not use 
this approach. Rather a two-layer model was used to integrate the equations all the way down 
to the wall. Sykes et al . (1986) and Oh and Schetz (1990) avoided this problem entirely by 
using slip conditions at the wall. Along symmetry planes, zero normal velocity and zero normal 
gradients for other variables are usually prescribed. Additionally, Reynolds shear stresses with a 
component in the normal direction will also be zero. Along the free stream, known variable values 
or zero surface stresses are prescribed. A major uncertainty exists as to the proper boundary 
condition at the jet exit plane. Experimental data by Andreopoulus (1983) had indicated that 
exit conditions are highly modified by the cross flow, especially at low R. Demuren (1983) 
found that by specifying constant total pressure at the jet exit, axial velocity profiles, similar to 
those observed experimentally could be simulated. However, in-plane velocity profiles had to be 
prescribed empirically. Kim and Benson (1992) overcame the uncertainty by placing the inflow 
boundary one diameter into the pipe from which the jet flow originates, and fully developed 
pipe flow conditions were prescribed at the recessed boundary. A different treatment would be 
required if the jet exited through an orifice rather than a pipe. 

Equations (26) - (41) form closed sets which can be solved by a finite-difference, finite 
volume or finite-element method to yield the mean flow and turbulence fields. By far the 
most popular approach is a combination of finite difference and finite volume methods. These 
are different manifestations and extensions of numerical techniques originally proposed by 
Chorin (1968), and Patankar and Spalding (1972). Notable exceptions are the finite-element 
computations of Baker et al . (1987) and Oh and Schetz (1990). 

Current numerical models can, in principle, be used to predict most flows of jets in cross 
flow which occur in practice. Both laminar and turbulent flows can be computed, so long as 
the flow is not dominated by rapid distortion or coherent structures. Instabilities develop in jet 
flows for Reynolds numbers greater than O (10), so that in most practical situations the flow will 
be turbulent or transitional. The former can be handled by current models with proper choice 
of the turbulence model. In the latter, large scale coherent structures may play a dominant role 
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[Andreopoulos (1985)], and methods based on Reynolds averaging may be inadequate. Direct 
calculation of the unsteady flow will generally be required, but this cannot be done for any 
practical Reynolds number, so a large eddy simulation appears to be a viable option. Such 
approaches are reviewed in Galperin and Orszag (1993), but much research work is still required 
before they can become reliable predictive tools in applications of interest here. In addition, 
current turbulence models, cannot predict higher-order statistics such as the Reynolds stresses 
[Speziale (1994)] in flows with rapid distortion (ratio of turbulent-to-mean-flow time scales 
greater than 50), or in flows with high compressibility (free shear flows with turbulent Mach 
number of order 1). Of course, the question of the importance of such higher-order statistics in 
the present flows has been raised. Some studies indicate that they are less important at high R, 
where pressure effects dominate, than at low to medium R. 

Earlier model computations of the fully elliptic type such as by Patankar et al. (1977) used 
grids (15x15x10 in the x, y, z directions) which are too coarse for the results to be considered 
reliable. Although correct trajectories were predicted, velocity fields deviated qualitatively and 
quantitatively from experimental observations. Finer grid computations have been reported, but 
grid independence could not be demonstrated conclusively in any of these. Multigrid methods 
[Demuren (1992)] allow systematic studies of grid dependence since on each finer grid level 
there are twice as many points in each direction than on the coarser one, i.e. eight times as many 
total grid nodes in three dimensions. Although up to 2.4 million nodes (i.e., 256x96x96 on the 
5th grid-level) were used, Claus and Vanka (1992) could still not demonstrate grid independence 
of the computed velocity and turbulence fields. Figure 4 shows comparisons of vertical profiles 
of the streamwise velocity component computed on the three finest grids. It is obvious that grid- 
convergence has not been achieved, especially in the near-field. Similar comparisons are shown 
for the turbulence level in Fig. 5. Far-field results appear closer to grid-convergence. With 
this type of grid refinement, better estimate of the results can be obtained by using Richardson 
extrapolation techniques [Demuren and Wilson (1994)]. The question of grid resolution cannot 
be completely separated from that of the formal order of accuracy of the numerical scheme. Most 
studies of the Patankar and Spalding (1972) type have used the hybrid upwind/central difference 
scheme to approximate convection terms. This is monotonic and conservative, but it is known to 
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be highly diffusive. [Demuren (1985b)]. Studies in which higher-order differences, such as the 
quadratic-upstream-weighted (QUICK) scheme were utilized e.g. Barata et al. (1991)] showed 
that similar results, as with lower-order schemes could be obtained on coarser grids. However, 
higher-order schemes tend to suffer from lack of boundedness in regions with high gradients. 

The uncertainty in the specification of boundary conditions for the jet hole exit has been 
discussed. Figure 6 shows contours of the jet velocity, static pressure and total pressure at the 
exit plane computed by Kim and Benson (1992). None of these is uniform, which is a compelling 
reason for including the jet pipe hole in the calculation domain. Other uncertainties involve the 
specification of inflow and near-wall conditions for turbulent quantities. 

In many computational studies, the inadequacy of the turbulence model has been blamed for 
the lack of agreement between computed results and experimental data. Demuren (1992) tried to 
isolate the effect of the turbulence model by performing computations on the same grid with the 
k-e model and the Reynolds stress (LRR-DH) model. The results are compared to experimental 
data of Atkinson et al. (1982) for opposed jets in cross flow in Figs. 7 and 8, for the mean flow 
and Reynolds stress fields, respectively. For the mean flow, there is little to choose between both 
model predictions, but the Reynolds stress model clearly gives better predictions of Reynolds 
stress profiles. From these results, it may be concluded that the mean flow was not strongly 
influenced by the turbulence field. Thus, if the interest is solely in the mean flow field, the k-e 
turbulence model, or an even simpler model, would be adequate. But, if the turbulence field is 
required, e.g., to predict mixing, then a Reynolds stress model would give much better results, 
but at additional computational cost. The multigrid technique enables the additional cost to be 
minimized by ensuring nearly grid-independent convergence rate. Thus, for a 3-level multigrid 
scheme, with (42x34x82) points on the finest grid, convergence of the Reynolds stress model 
computations could be obtained in less than 100 fine grid iterations. 

Concluding Remarks 

In the study of turbulent jets in cross flow, empirical models offer quick and simple methods 
for obtaining first-order estimates and a qualitative picture of the jet trajectory, its extent, and 
decay rates of its axial velocity and temperature. The main requirement for reasonable predictions 
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is the use of correlation equations or curves derived from experimental data bases with similar 
characteristics as the problem of interest. 

Integral models contain in simplified forms mathematical representations of the basic con- 
servation laws, and can thus be applied much more widely than empirical models. Several 
physical phenomena which occur in the flow are modeled with relations which are more or less 
empirical. Combinations of these have been used successfully in integral models, so long as 
they are properly calibrated. One criticism of integral models is that they provide little insight 
into flow physics, since the same effect could be achieved in several different ways. All integral 
models are computationally cheap to use. The basic models are conceptually simple, but more 
sophisticated models have been devised which enable more complex jet cross flow interactions 
to be analyzed. 

Perturbation models do not require much empirical input, but they are mostly restricted to 
the near-field or far-field where small parameters required for expansions can be defined. They 
enable order-of-magnitude studies of the effects of various parameters, and are thus useful tools 
for the investigation of flow physics. Beyond these, they have limited practical utility. 

Numerical models offer the best choice as practical predictive tools over a wide range of jets 
in cross flow applications. They require the least assumptions and empirical input. They are, 
however, the most computationally intensive. Quite complex jet-jet, jet-cross flow interactions 
can be analyzed. Depending on specific requirements, the choice of turbulence model may or 
may not be important. Although, several quite complex flow situations have been computed with 
a measure of quantitative accuracy, some questions remain as to the effects of grid resolution, 
turbulence model and boundary conditions on overall accuracy of computed results. Reliability 
and computational accuracy are expected to improve with further developments in numerical 
techniques and turbulence models. These are clearly the models of choice for the computation 
of practical jets in cross flow situations. Computer codes are available commercially for this 
purpose. 
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FIGURE CAPTIONS 


Figure 1 : Configuration of a jet in cross flow 

Figure 2 : Prediction of jet trajectories and concentration decay; symbols — experimental 

data (from Fan, 1967), curves — calculations (from Abraham, 1971). 

Figure 3 : Nomenclature for perturbation model description (from Coelho and Hunt, 1989); 

(a) 2D vortex sheet model, (b) 3D vortex sheet model, (c) entraining 3D vortex 
sheet model. 

Figure 4 ; Grid dependency test: comparison between computations from Claus and Vanka 
(1992) and experimental data from Khan et al. (1982) of stream wise velocity 
in the center plane; (a) x/D = 4, (b) x/D = 8. 

Figure 5 : Grid dependency test: comparison between computations from Claus and Vanka 

(1992) and experimental data from Crabb et al . (1981) of turbulence intensity 
in the center plane; (a) x/D = 2, (b) y/D = 1.35. 

Figure 6 : Contours plots of the flow field at the jet exit plane (from Kim and Benson, 

1992); (a) axial velocity (Uj/Oj), (b) static pressure (P / ^pUj), (c) total pressure 

[(p+ipufyipufi. 

Figure 7 : Effect of turbulence model: comparison between computations from Demuren 

(1992) and experimental data from Atkinson et id. (1982) of streamwise velocity 
in the center plane; (a) flow configuration - opposed jets in cross flow, (b) x/D 
= 8, (c) x/D = 12. 

Figure 8 : Effect of turbulence model: comparison between computations from Demuren 
(1992) and experimental data from Atkinson et al. (1982) of Reynolds stresses 
in the center plane; (a) x/D = 8, (b) x/D = 12. 
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